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Matter described by Quantum Chromodynamics (QCD), the theory of strong 
interactions, may undergo phase transitions when its temperature and the 
chemical potentials are varied. QCD at finite temperature is studied in the 
laboratory by colliding heavy-ions at varying beam energies. We present a 
test of QCD in the non-perturbative domain through a comparison of thermo- 
dynamic fluctuations predicted in lattice computations with the experimental 
data of baryon number distributions in high-energy heavy-ion collisions. This 
study provides evidence for thermalization in these collisions, and allows us to 
find the crossover temperature between normal nuclear matter and a decon- 
fined phase called the quark gluon plasma. This value allows us to set a scale 
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for the phase diagram of QCD. 

QCD is the theory of strong interactions — one of the four fundamental interactions oc- 
curring in nature, and an essential part of the standard model of particle physics. It describes 
interactions between quarks and gluons, which are the ultimate constituents of the majority of 
the visible mass of the universe (1, 2). In the short-distance regime where the momentum ex- 
change between quarks and gluons is large, the strong coupling constant becomes small through 
the mechanism of asymptotic freedom. In this perturbative region QCD is very successful in ex- 
plaining various processes observed in experiments involving electron-positron, proton-proton 
and proton- antipro ton collisions (3). In the non-perturbative regime tests of the theory were 
related to the computation of hadron properties (4). In other regimes of long-distance non- 
perturbative physics, the theory is yet to be tested. Here, we test the thermodynamics of bulk 
strongly interacting matter. 

Experimental tests of non-perturbative QCD in the bulk can be carried out by colliding 
heavy-ions (like U, Pb, Au, and Cu) at different center of mass energies, y/s NN (5-8). Several 
experimental programs have been launched or are in the planning stage at facilities such as the 
Large Hadron Collider (LHC), the Relativistic Heavy Ion Collider (RHIC), the Super Proton 
Synchrotron (SPS), the Facility for Anti-proton and Ion Research (FAIR) and the Nuclotron- 
based Ion Collider f Acility (NICA), where the essential features of the QCD phase diagram can 
be studied. 

In QCD there are conserved quantities like the net-baryon number, B, the net-electric 
charge, Q, and the net-strangeness, S. The term net means the algebraic sum of the quan- 
tum numbers, where those of anti-particles are the negatives of the corresponding particles. As 
a result the thermodynamics of the bulk can be characterized by the corresponding chemical 
potentials (energy needed to add/remove one unit of the conserved quantity to/from the system) 
li B , /jl q , and /ji s in addition to the temperature, T, conjugate to the conserved energy of a bulk 
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system. In experimental studies of particle ratios measured in heavy-ion collisions it is observed 
that the relevant values of \i Q and fj, s are small compared to ji B . For example, in Au ion colli- 
sions within rapidity range of ± 0. 1 unit at y/s NN = 200 GeV (with impact parameter less than 
3 fm) one finds that /i B = 22 ± 4.5 MeV, while /i s = 3.9 ± 2.6 MeV and fi Q is still smaller (9). 

The lattice formulation of QCD is a non-perturbative approach from first principles for 
obtaining the predictions of QCD. Space-time is replaced by a lattice; quarks occupy the sites, 
and gluons occupy the links between the sites. The lattice spacing, a, is the inverse of the cutoff 
required to regulate any interacting quantum field theory. The theory is solved numerically at 
several values of a. The extrapolation to the continuum (a = 0) can then be made through the 
renormalization group equations. In QCD there is a conventional temperature, T c , which is an 
intrinsic scale of bulk hadronic matter. We follow the definition that it is the temperature at the 
peak of a susceptibility related to the confinement-deconfinement order parameter (called the 
Polyakov loop susceptibility, Xl) at /i s = (10-13). Lattice QCD computations show that this 
peak is finite, which corresponds to a cross-over (14, 15). The temperature at which Xl peaks, 
of course, changes with /i s . However, once T c is known, such shifts as a function of /x s can 
be quantified. This is similar to saying that the Celsius scale of temperature is defined by the 
boiling point of water at normal pressure, P, and that the boiling point changes with P. 

One of the most basic questions to ask about bulk hadronic matter is the value of T c . This 
can be represented as a link in a "circle of reasoning" that encompasses all the regimes of 
non-perturbative QCD (Fig.l). So far the strategy to find T c has been indirect: first lattice QCD 
computations are performed at both T = and T > in order to determine a ratio T c /m, where 
m is a typical hadronic scale (step (b) of Fig.l). Then one replaces the scale m, determined on 
the lattice, by an experimental measurement (step (a)). The temperature at each v /^w extracted 
from models of particle yields (18, 19) is step (d) of the circle of reasoning. From such models 
one finds that the fireball of bulk nuclear matter created in heavy-ion collisions, which is initially 
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Figure 1 : Illustration of the chain of reasoning for testing QCD in the non-perturbative domains 
of the strong interactions and obtaining the scale, T c , of the QCD phase diagram. 

out of equilibrium, evolves to a state of thermal equilibrium at chemical freeze-out. The models 
do no give T c ; however, they allow the extraction of T and ji B at freezeout. In this paper we show 
that predictions of lattice computations of finite temperature QCD (16), taken in conjunction 
with determinations of T c in step (b) (10-13) agree well with experimental measurements on 
bulk hadronic matter (17). This allows us to invert the reasoning and extract T c directly from the 
experimental measurements in heavy-ion collisions (step (c)). The agreement of the temperature 
from steps (c) and (d) along with the agreement of T c extracted from steps (a) and (b) with that 
from (c) show the complete compatibility of a single theory of hadron properties and of bulk 
QCD matter, i.e., of all non-perturbative regimes of the strong interactions. This approach may 
present a new domain of tests of the standard model of particle physics. 

The conjectured phase diagram of QCD: In the current conjectures for the parts of the 
phase diagram that is accessible with heavy-ion collisions (20) (Fig. 2), calculations within sim- 
plified models which mimic QCD show that at large fi B there is a first order hadron-QGP phase 
transition. This phase boundary is expected to end in a critical point at finite [i B , as lattice com- 
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Figure 2: Current conjectures for the QCD phase diagram. The phase boundary (solid line) be- 
tween the normal low-temperature hadronic phase of bulk QCD matter and the high-temperature 
partonic phase is a line of first order phase transitions which begins at large fx B and small T and 
curves towards smaller ji B and larger T. This line ends at the QCD critical point whose probable 
position, derived from lattice computations, is marked by a square. At even smaller fj, B there are 
no phase transitions, only a line of cross-overs (shown by a dashed line). The red-yellow dot- 
ted line corresponds to the chemical freeze-out line from the evolution of the bulk QCD matter 
produced in high energy heavy-ion collisions. The solid point at T = and [i B = 938 MeV 
represents nuclear matter in the ground state. At large /i B and low T is the color superconductor 
phase (CSC) (26). 
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putations (10-13) agree with general symmetry arguments (21), which indicate that at /j, B = 
there is neither a first-order nor a second-order phase transition but only a cross-over at T c . 
The determination of T c sets the scale of the QCD phase diagram. Current best estimates of 
the position of the critical point (22) are reflected in the position indicated in Fig. 2. The ex- 
perimental focus currently is on an attempt to locate the critical point and the line of phase 
coexistence (23, 24). 

By changing y/s NN one traces out a line of chemical freeze-out in the phase diagram, as 
shown in Fig. 2. This line is parameterized through a hadron resonance gas model (18, 19). 
Because this work focuses on making a connection between QCD thermodynamic calculations 
and observables measured in experimental facilities, we also show in Fig. 2 the range of /i s 
values covered by various experiments as one traverses the chemical freeze-out line by changing 
^s NN . The solid point around fi B = 938 MeV is the location of ordinary nuclear matter (25), 
the best characterized point on the phase diagram. 

Comparison of experimental measurements with lattice QCD predictions: Lattice QCD 
computations leave open the question of a scale and yield dimensionless predictions, for exam- 
ple for P/T 4 as a function of T/T c and /jl b /T. Here we discuss the non-linear susceptibilities 
(NLS) of baryons, x^b \ of order n (27). These are the Taylor coefficients in the expansion of P 
with respect to fx B at fixed T in the usual dimensionless form 



Lattice measurements of the series expansion of the NLS in powers of fi B /T are resummed 
using Pade approximants in order to give predictions for the above quantities (16). They are 
of interest because they are related to cumulants of the fluctuations of the baryon number in 
thermal and chemical equilibrium in a grand canonical ensemble. 




(1) 
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The n-th cumulant of such fluctuations, [£>"], is given by 

[B n } = VT S T n^ x (n) ^ ^ ? (2) 

where V is the volume of the observed part of the fireball. Because observed hadrons are in 
thermal and chemical equilibrium at the freeze-out, this relation should hold for cumulants 
of the observed event-by-event distribution of net-baryon number in heavy-ion collisions. The 
cumulants are often reported as the variance, a 2 = [B 2 ], the skewness, S = [B 3 ]/[B 2 ] 3 ^ 2 and the 
Kurtosis,/? = [_B 4 ]/[i? 2 ] 2 . It is clear from these definitions that the ^-dependence in Eq. [2]gives 
the correct volume scaling expected from the central limit theorem. This leads to the classic 
extraction of the susceptibility from fluctuations in the grand canonical ensemble (28, 29). 

There is one remaining subtlety in comparing lattice computations with experimental data. 
Most experiments are designed to measure event-by-event net-protons. The data discussed in 
the current work is from the STAR experiment at RHIC (17), which identifies protons and anti- 
protons by measuring the specific ionization energy loss of these particles in the gas of a Time 
Projection Chamber. These measurements miss neutrons, the other dominant part of the baryon 
distribution. This may impose limitations on our measurement of fluctuations. However the 
effect of isospin fluctuations on the shape of the net-baryon distributions is small (30). Hence 
we proceed under the assumption that the shape of the net-proton distributions reflects the net- 
baryon distributions up to distortions smaller than the estimated errors in measurements of the 
cumulants. 

We are unable to exploit Eq. |2]directly in heavy-ion experiments because the volume, V, is 
hard to determine precisely experimentally. However, the ratios 

[B 3 ] T X f 



(mi) : Sa 



[R 4 1 

(m 2 ) : kg 



2 _ [£ 4 ] _ T 2 X L 4) 
" IB 2 ] xi 2) 
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K) : -j = W] = -^p 0) 

do not contain the volume and therefore provide a direct and convenient comparison of ex- 
periment and theory (31). The above equations are written in a form that emphasizes this 
connection — the left hand side can be measured in an experiment whereas the right hand side 
can be predicted by lattice QCD. We use the notation m 1 2j 3 generically to refer to either side. 

We now proceed to discuss the comparison of m\ and m 2 from experiment and theory 
(Fig. 3). The experimental measurements (17) were made using the number of protons (p) and 
anti-protons (p) produced in the collision of Au ions around 90° to the beam axis with the impact 
parameter of the collisions being less than 3 fm (32). p and p are in the range of 400 MeV/c 
to 800 MeV/c. This choice of momentum range is designed to obtain the purest sample of p 
and p. A large fraction of p and p is contained in this kinematic range. The effect of finite 
reconstruction efficiency of p and p has been shown to be negligible (17). The experimental 
values of Sa and na 2 are shown as a function of y/s NN . 

The lattice calculations (16) were carried out using two flavors of staggered quarks in QCD. 
The lattice cutoff 1/a ~ 960 — 1000 MeV and the bare quark mass were tuned to give a pion 
mass of about 230 MeV (33). These computations were performed at p, B = and the Taylor 
series coefficients of P/T 4 were used to extrapolate mi and m 2 to the freeze-out conditions 
using appropriate order Pade approximants to resum the series expansions. Since lattice results 
are obtained in terms of T/T c and fi B /T (see Eq. Q]), their extrapolation to the freeze-out 
conditions required the input of T c . The lattice values were obtained using T c = 175 MeV, 
compatible with indirect determinations of T c (10-13). 

On the upper scales of Fig. 3 we also show the p, B and T values at chemical freeze-out that 
correspond to the various y/s NN . For this we used the functional relationship between these 
values from the hadron resonance gas model using the yields of hadrons discussed in (18, 19). 



8 



(A) 



CO 



(B) 



79 



140 



720 420 



T (MeV) 
160 

u (MeV) 
210 54 



165 



166 



20 10 



1.5 



1 

0.5 

2 
1 

* 
-1 
-2 



■ 1 ' — ' — ' — 1- >- 1 r 

- m 7 • l 

L "'$ 
: i > 


1 1 . 1 1 1— i— | 1 1 " 

Exp. Data Lattice QCD - 
HRG \ 

T - 

" 

1 — | 1 1—: 


i m 2 




41 


¥ ¥ • 

: 

{ ; 

i ■ ■ — 



CM 



CO 



CM 



CM 



4 56 10 20 30 100 200 

Vsnn (GeV) 



Figure 3: Comparison of lattice QCD and experimental data for mi (A) and m 2 (B). Experi- 
mentally measured ratios of cumulants of net-proton distributions, m x = Sa and m 2 = na 2 are 
shown as a function of y/s NN for impact parameter values of less than 3 fm for Au+Au colli- 
sions at RHIC (17). Also plotted on the top scale are the chemical freeze-out values of fx B and 
T corresponding to y/s NN as obtained from a hadron resonance gas model, which assumes the 
system to be in chemical and thermal equilibrium at freeze-out (18, 19). The prediction of such 
a model for mi (34) is shown as the dashed red line. The lattice predictions for these quantities 
are drawn from a computation with lattice cutoff of 1/a ~ 960 — 1000 MeV and converted to 
the dimensional scale of T and [i using T c = 175 MeV. 
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The model predictions for mi (34) are also shown. The hadron resonance gas model predictions 
can be reproduced if baryon and anti-baryon numbers are independently Poisson distributed. 
Having established a connection between y/s NN and (T, fx B ) we compare the experimental data 
on fluctuations to those predicted from lattice QCD. Excellent agreement is seen between lattice 
QCD predictions and experimental measurements for all three beam energies. This marks the 
first successful direct test of QCD against experimental data in the non-perturbative context of 
bulk hadronic matter. The agreement with the data is yet another non-trivial indication that the 
fireball produced in heavy-ion collisions is in thermal and chemical equilibrium at chemical 
freeze-out. 

Setting the scale of bulk QCD: Lattice QCD results for mi ^3 are obtained for dimension- 
less arguments T/T c and ji B /T as shown in Eq. 2. For a given value of ^/s NN , the experimental 
observations are realized at the corresponding chemical freeze-out, characterized by T and fx B . 
Thus a comparison of the two requires a choice of the scale, T c . By varying this scale to obtain 
the best fit between the QCD predictions and experimental measurements, we determine T c for 
the first time through an observable connected to strongly interacting bulk matter. The results 
are, of course, subject to all the caveats expressed in the previous section. The observable that 
we choose for comparison is m 3 . The lattice computation of this quantity has the smallest 
systematic uncertainties among the three explored here, and thus is the best quantity to use to 
constrain T c . 

Figure 4A shows the comparison of m 3 between experimental results from Au ion collisions 
and lattice QCD predictions. This is an extension of Fig. 3 which shows a comparison with mi 
and m 2 . In this analysis, the results of m l5 m 2 and m 3 are consistent as required in Eq. |3] 
The new information here is that we show lattice predictions obtained with different values of 
T c . The errors on the experimental data points are statistical (lines) and systematic (brackets) 
errors (17). The errors bars on the lattice predictions are statistical errors with a cutoff of 
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Figure 4: Comparison of m 3 from experiment and lattice predictions, and the extraction of T c . 
(A): ku/S of net-proton distribution measured in collisions of Au ions at varying y/s NN and 
with an impact parameter of less than 3 fm. This is compared to lattice QCD predictions with 
cutoff 1/a ~ 960 — 1000 MeV for the corresponding ratio of susceptibilities extrapolated to the 
freeze-out conditions using different values of T c . The lattice results at each y/s NN are slightly 
shifted for clarity in presentation. (B): The comparison of experimental data and lattice QCD 
predictions, shown through x 2 as a function of T c using the definition given in Eq. |U This 
yields the estimate of T c and its errors as discussed in the text. 
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1/a ~ 960 — 1000 MeV. The lattice spacing effects and the effect of tuning the bare quark mass 
are the main sources of remaining uncertainties in the predictions. These are not parameterized 
as systematic uncertainties. However, it is known that their effect is small at the two highest 
values of y/s NN (16). 

In order to arrive at a quantitative estimate of the scale parameter T c we perform a standard 
statistical analysis. For each value of T c we compute, 

\ — V" [ m 3 P {\/ s nn) — (V-Sjviv, T c )] 2 ^ 

V*im Error ex P t + Error^ CD 

where the errors in the experimental and lattice QCD quantities are obtained as explained above. 
The lattice predictions are calculated for the grid of T c values (Fig. 4). The minimum of x 2 > 
corresponding to the most probable value of the parameter being estimated, occurs at T c = 
175 MeV. The standard errors on the parameter are the values of T c for which x 2 exceeds the 
minimum value by unity. It is clear from Fig.4(B) that this is bounded by +5 and —10 MeV. A 
piece-wise linear interpolation between the grid points yields the more reliable error estimate, 
+1 and —7 MeV. By comparing different interpolation schemes we find that the error estimate 
is stable. As a result we conclude that 

T c = 175±\ MeV. (5) 

The error estimates include systematic and statistical errors on experimental data but only sta- 
tistical errors on the lattice QCD computations. 

The result in Eq. [5] is compatible with current indirect estimates of T c which come from 
setting the scale of thermal lattice QCD computations via hadronic observables. Furthermore, 
this gives a scale for temperatures which is compatible with the resonance gas model, as shown 
in Fig. 3. As we discussed in the introduction, this closes a circle of inferences which shows that 
phenomena obtained in heavy-ion collisions are fully compatible with hadron phenomenology, 
and provides a first check in bulk hot and dense matter for the standard model of particle physics. 
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Conclusions and Outlook: We have performed a direct comparison between experimental 
data from high energy heavy-ion collisions on net-proton number distributions and lattice QCD 
calculations of net-baryon number susceptibilities. The agreement between experimental data, 
lattice calculations and a hadron resonance gas model indicates that the system produced in 
heavy-ion collisions attained thermalization during its evolution. The comparison further en- 
ables us to set the scale for non-perturbative high temperature lattice QCD by determining the 
critical temperature for the QCD phase transition to be 175 MeV. 

This work reveals the rich possibilities that exist for a comparative study between theory 
and experiment of QCD thermodynamics and phase structure. In particular, the current work 
can be extended to the search for a critical point. In a thermal system, the correlation length (£) 
diverges at the critical point. £ is related to various moments of the distributions of conserved 
quantities such as net-baryons, net-charge, and net-strangeness. Finite size and dynamical ef- 
fects in heavy-ion collisions put constraints on the values of £ (35). The lattice calculations 
discussed here and several QCD-based models have shown that moments of net-baryon distri- 
butions are related to baryon number susceptibilities and that the ratio of cumulants m 2 = no 2 , 
which is related to the ratio of fourth order to second order susceptibilities, shows a large devi- 
ation from unity near the critical point. Experimentally, na 2 can be measured as a function of 
^/s NN (or T and /j, b ) in heavy-ion collisions. A non-monotonic variation of no 2 as a function 
of y/s NN would indicate that the system has evolved in the vicinity of the critical point and thus 
could be taken as evidence for the existence of a critical point in the QCD phase diagram. 
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